arXiv:1506.08985vl [hep-ex] 30Jun2015 


Systematic study of the experimental measurements on J/ip cross section and 
kinematic distribution in p + p collisions at different energies 

Wangmei Zha, 1 ’ 2 Bingcliu Huang, 2 Rongrong Ma, 2 Lijuan Ruan, 2 Zebo 
Tang, 1 ’^] Zhangbu Xu, 2 Chi Yang, 1 ’ 2 Qian Yang, 1,2 and Shuai Yang 1,2 

1 University of Science and Technology of China , Hefei, China 
2 Brookhaven National Laboratory, New York, USA 
(Dated: July 1, 2015) 

The world experimental data on cross section and kinematic distribution in p+p and p+A collisions 
at yjs = 6.8 - 7000 GeV are examined in systematic way. The yjs dependence of the inclusive 
cross section, rapidity and transverse momentum distributions are studied phenomenologically. We 
explore empirical formulas to obtain the total cross section, rapidity and transverse momentum 
(Pt) distribution. This is crucial for the interpretation of A+A J/ip results at RHIC when the 
p + p reference data are not available. In addition, the cross section at mid-rapidity and transverse 
momentum distributions in p + p collisions at y/s = 39 and 62.4 GeV are evaluated. 


I. INTRODUCTION 

Lattice QCD predicts that, under conditions of ex¬ 
tremely high temperatures and energy densities, a phase 
transition or crossover from hadronic matter to a new 
form of matter, known as Quark Gluon Plasma (QGP) 
P], will occur. The Relativistic Heavy Ion Collider 
(RHIC) was built to search for the QGP and to study its 
properties in laboratory through high-energy heavy-ion 
collisions [2). Many observables have been proposed to 
probe the QGP created in heavy-ion collisions. Among 
them, the J/ip suppression caused by the color-charge 
screening in QGP is one of most important signatures 

El- 

Over the past twenty years, J/ip production in hot 
and dense medium has been a topic attracting growing 
interest. Suppression of J/ip production has been ob¬ 
served in various experimental measurements @HZj. A 
similar suppression pattern and magnitude of J/ip was 
observed at SPS and RHIC despite of huge collision en¬ 
ergy difference. Furthermore, the J/if is suppressed more 
in forward rapidity than that in midrapidity at RHIC 
200 GeV Au+Au collisions JHj . These experimental ob¬ 
servations suggest that, in addition to color screening, 
there exist other effects contributing to the modification 
of J/ip production. Cold nuclear matter (CNM) effects, 
the combined contribution of finite J/ip formation time 
and finite space-time extent of QGP and recombination 
from uncorrelated c and c in the medium may account for 
these contributions [S]- Among these contributions, the 
regeneration of J/ip from the recombination of cc plays 
an important role to explain the similar suppressions at 
SPS and RHIC. As the collision energy increases, the re¬ 
generation of J/ip from the larger charm quark density 
would also increase which partly compensates for the ad¬ 
ditional suppression from color-screening. The regenera¬ 
tion also expects a stronger suppression at forward rapid¬ 
ity at RHIC where the charm quark density is lower than 
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that at midrapidity. At LHC, the J/ip is less suppressed 
in both mid-rapidity and forward rapidity than that at 
RHIC TO) HI], which may indicate that the regenera¬ 
tion contribution is dominant in the J/ip production at 
LHC energies. Measurements of J/ip in different collision 
energies at the Solenoidal Tracker at RHIC (STAR) can 
give us indications on the balance of these mechanisms 
for J/ip production and medium properties. 

To qualify the medium effects on the modification of 
J/ip production, the knowledge of J/ip cross section and 
kinematics in p+p collision is crucial to offer a reference. 
During RHIC year 2010, STAR has collected abundant 
events of Au+Au collisions at Jsnn = 39 and 62.4 GeV, 
while the reference data in p + p collisions is not in the 
schedule of RHIC run plan. As what we did in ref fl2l . 
we study the world-wide data to obtain the J/ip reference 
at these collision energies. 

In this letter,we report an interpolation of the pp- 
integrated and differential inclusive J/ip cross section in 
p+p collisions at mid-rapidity to yjs = 39 and 62.4 GeV. 
We establish a strategy to estimate the inclusive J/ip 
cross section and kinematics at certain energy points, 
which makes the calculation of the J/ip nuclear modi¬ 
fication factors for any colliding system and energy at 
RHIC possible. The extrapolation is done in three steps. 
The first step is an energy interpolation of the existing 
total J/ip cross section measurements. The second step 
is the description of the energy evolution of the rapidity 
distribution. The last step is the evaluation of the energy 
evolution of the transverse momentum distribution. 


II. AVAILABLE EXPERIMENTAL RESULTS 
TREATMENT 

The measurements of J /ip hadroproduction have been 
performed for about forty years. In such a long period, 
different experimental techniques have been utilized and 
different input information was available at the time of 
the measurements. Therefore, comparison of different 
experimental results on an equal footing needs an up- 
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date of the published values on several common assump¬ 
tions and aspects. For example, the branching ratio of 
J/ip —>• e + e _ (or p + p~) have changed with time; the as¬ 
sumed functional forms for the xf and pt shapes, which 
can be used to infer the total J/ip production, are differ¬ 
ent in different measurements; and the treatment of the 
nuclear effects are not homogeneous. In this section, we 
update all the results with the current best knowledge of 
branching ratios, kinematics and nuclear effects. 

The cross section for J/ip on a nuclear target is often 
characterized by a power law: 

X r- (1) 

where is the J/ip proton-nucleon cross section and 

<j J Ja is the corresponding proton-nucleus cross section 
for a target of atomic mass number A. 



FIG. 1. (color online) Measurements of a as a function of 
xf by various experiments in different collision energies. The 
solid curve represents the parametrization discussed in the 
text. 


The dependence of a on xf measured by NA3 m 
NA50 m E772 US], E886 [IS] and HERA-B [T7j is 
shown in Fig. [l] where Xp is defined as Xp = 2 p z /yfs ( p z 
is longitudinal momentum, along the beam direction.). 
No significant energy dependence of a as a function of 
xf is observed within uncertainties, thus we assume it 
is independent of the cms-energy (y^). The results of 
J/il> a at Xf > 0 can be represented for convenience by 
the simple parametrization shown as solid line in Fig. |Tj 
cx(xf) = 0.9503e -in2 (i 384e) . If an experiment has 

published cross sections of J/ip in proton nucleus col¬ 
lisions, Eq. [T] and the solid curve in Fig. [l] is applied 
to obtain the corresponding J/'ip cross section in proton 
nucleon collisions. Some of the experimental measure¬ 
ments are only quoted for a limited phase-space. To ob¬ 
tain the total cross sections, the functional forms of xp 
and pt spectrum shapes utilized for extrapolation are: 
da/dxF = a x e ~ ln2( -~^~ E ,and da/dpp = a x w 

respectively, where a, b, and c are free parameters. As 


illustrated in Fig. [2] these two functional forms describe 
the Xf and pp spectra very well. All the measure¬ 
ments are updated with the latest branching fractions 
(5.961 ± 0.032% for J/ip -> p+ + p~ , 5.971 ± 0.032% for 
J/ip — e + + e~) [TS]. The treated results on J/ip cross 
sections are listed in Tab.[I]They show a good overall con¬ 
sistency, even though some of them contradict with each 
other. For example, the two measurements (E331 and 
E444) at 20.6 GeV deviate from each other by roughly 
2cr. The E705 measurement at 23.8 GeV is higher than 
the UA6 one at 24.3 GeV by more than 2 ct. 

III. RESULTS 

The energy evolution of the total inclusive J/ip produc¬ 
tion cross section in proton induced interactions is shown 
in Fig. [3] The first approach is to use the predicted shape 
in the Colour Evaporation Model at Next to Leading Or¬ 
der (NLO) [40] to describe the energy dependence of 
J/ip cross section. The central CT10 parton density set 
[TT1 and {m, pf/tu, pp/m} = {1.27(GeV), 2.10,1.60} set 
is utilized in the predicted shape, where m is the charm 
quark mass, pp is the factorization scale, pp is the renor¬ 
malization scale. The fit is defined such that the nor¬ 
malization of the NLO CEM calculation is left as a free 
parameter (a): a = a x ijcem■ The second approach 
is to use a functional form to describe the cross section 
energy evolution: f{y/s) = a x y!/ lax x e«~'+=, where 
Vmax = ln( ~ys )) a, b, c and d are free parameters. As 

J / ip 

shown in Fig. [ 3 ] both approaches can describe the en¬ 
ergy evolution trend of J/ip cross section. The y 2 /NDF 
for CEM and functional fit are 90.5/22 and 76.7/20, re¬ 
spectively. The large y 2 mainly comes from three exper¬ 
imental points which contradict with the common trend 
(E331 and E444 measurements at 20.6 GeV, E705 mea¬ 
surement at 23.8 GeV). If we exclude these three data 
points and refit the results, the y 2 /NDF for CEM and 
functional fit are 41.1/19 and 16.7/16, respectively. The 
values extrapolated (without the three bad experimental 
points) for the J/ip cross sections at yfs = 39 and 62.4 
GeV, utilizing the functional form and the NLO CEM 
based fit are listed in Table mi 

The knowledge of the rapidity dependence of J/ip pro¬ 
duction at different cms-energies is crucial to obtain a ref¬ 
erence for the measurements at mid-rapidity from RHIC. 
Based on a universal energy scaling behavior in the ra¬ 
pidity distribution obtained at different cms-energies, we 
explore approaches to the extrapolation of the rapidity 
distribution. As shown in Fig. [4] the y-differential cross 
sections at different cms-energies have been normalized 
by the total cross section, and the normalized values 
are plotted verse y/y m axi where y max has been previ¬ 
ously defined. Despite of huge cms-energy difference, the 
treated RHIC [3B] and LHC [381 [Ml H2j experimen¬ 
tal distributions fall into a universal trend, which allows 
us to perform global fits to all the experimental results 
with suitable functions. Two functional forms are cho- 
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FIG. 2. (color online) Distributions of Eda/dpr (left panel) and Edcr/dxF (right panel) in p + C collisions at y / s = 20.6 GeV 
measured by E331 collaboration |23| . The solid lines are fit curves with the functional forms described in the text. 


Experiment 

Reaction 

(GeV) 

<jj/$ (nb/nucleon) 

CERN-PS [19] 

p±A 

6.8 

0.732T0.13 

WA39 [20] 

P±P 

8.7 

2.35T1.18 

IHEP [21 

p±Be 

11.5 

21.63T5.64 

E331 [22] 

p±Be 

16.8 

85.15T21.30 

NA3 pJH 

p±Pt 

16.8 

95.0T17.0 

NA3 Q3] 

p±Pt 

19.4 

122.6T21 

NA3 Q3] 

P±P 

19.4 

120±22 

E331 [23j 

P±C 

20.6 

278±32.8 

E444 (1241 

P±C 

20.6 

176.5T23.3 

E705 [29] 

p±Li 

23.8 

271.51±29.84 

UA6 [2B] 

P±P 

24.3 

171.42T22.21 

E288 [27] 

p±Be 

27.4 

294.12T73.53 

E595 [28] 

p±Fe 

27.4 

264±56 

NA38/51 [30l|31] 

p±A 

29.1 

229.5±34.4 

NA50 [13] 

p±A 

29.1 

250.7T37.6 

E672/706 [32] 

pBe 

31.6 

343.07T75.12 

E771 [33] 

p±Si 

38.8 

359.1T34.2 

E789 [34] 

p±Au 

38.8 

415.04±100 

ISR [46] 

P±P 

52 

716T303 

PHENIX [3g] 

P±P 

200 

4000T938 

CDF [37] 

p+p 

1960 

22560±3384 

ALICE [38| 

P±P 

2760 

29912.6±5384.3 

ALICE [39] 

P±P 

7000 

54449.4±8494 


TABLE I. (color online) Updated total (aj/^,) production cross sections in proton-induced interactions. 


Fit 

cross section (nb/nucleon) 

yjs =39 GeV 

y/s =62.4 GeV 

NLO CEM 

425±20 

941±42 

function 

445±30 

995±66 

evaluated results 

445±30±20 

995±66±54 


TABLE II. Extrapolated values of the J/if) production cross 
section at = 39 and 62.4 GeV. The difference between 
CEM and function fit has been taken as the systematic un¬ 
certainties of the extrapolation. 


sen to do the fits: one is Gaussian function, the other 

one is 1 _( y /y -pe v u-v/in.d' , where a and b are 

free parameters. Both of them can describe the global 
distribution very wel \(x 2 /NDF = 10.1/27 for gaussian 
fit, x 2 /NDF = 11.2/27 for the other fit). With the 
extrapolated J/if) cross sections and rapidity distribu¬ 
tions, the predicted J/ip cross section times branching 
ratio at y/s = 39 and 62.4 GeV in mid-rapidity are 
Br(e + e~)da/dy| M<1 . 0 = 9.04 ± 0.69 and 17.74 ± 1.06 
nb, respectively. The uncertainties include statistical and 
systematic uncertainties. These values are highly consis¬ 
tent with the estimations from CEM model (8.7 ±4.5 nb 
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FIG. 3. (color online) Energy dependence of inclusive J/ip 
prodcution cross section. The open circle is the fit from CEM 
shape.The solid line is a function fit as discuss in the text. 



FIG. 4. (color online) Normalized J/i/> production cross 
section as a function of y/y ma x- Two function fits 
are shown: one is Gaussian function, the other one is 

i_^ y /y -p-e ( '!-»/>"«■. The difference between these 

two curves has been considered as systematic errors. 


for 39 GeV, 17.4 ± 8.0 for 62.4 GeV). 

The energy evolution of J/ij) transverse momentum 
distribution are also studied via available experimental 
measurements from y/s = 10 - 7000 GeV [TS] [531 [371 
1331133113711331133] . We use part of the world-wide fixed- 
target data (with only p, Be, Li, and C respectively) 
measured with incident protons. In this way, we avoid 
uncertainties due to ignoring any cold nuclear matter ef¬ 
fects on the J/'ip transverse momentum distributions. To 
compare the different experimental measurements at dif¬ 
ferent energies and rapidity domains, as shown in Fig. [5] 
the transverse momentum distributions are normalized 
by their pr-integrated cross sections and plotted verse the 
zt variable, which is defined as Zt = Pt/ < Pt >■ The 
treated distributions follow a universal trend despite of 
the different cms-energies and rapidity domains. We can 



FIG. 5. (color online) J/ij) zt distributions for available ex¬ 
perimental results from y/s = 10 to 7000 GeV. The solid line 
is a function fit as discussed in the text. 


describe the global distributions very well by the follow¬ 
ing function: = a x {1+ J z , )n 03], where 

a = 2 b 2 (n - 1), b = T(3/2)T(n - 3/2 )/T{n - 1), and 
n is the only free parameter. From the fit, we obtain 
n = 3.93 ± 0.03 with x 2 /NDF = 143.9/162. 



FIG. 6. (color online) J/ij) < pr > at mid-rapidity as a 
function of cms-energy from y 7 i = 10 to 7000 GeV. The solid 
line is a function fit as discussed in the text. 


With the universal shape and < pp > information at 
certain energy and rapidity domain (we focus on mid¬ 
rapidity), we can extrapolate the transverse momentum 
distribution at any cms-energy. Thus the next step is to 
evaluate the energy evolution of < pp >. The < pr > at 
mid-rapidity as a function of cms-energy from world-wide 
experiments US 121 EZl E21 ESI EH EH mMj is shown in 
Fig-© Also, only part of the world-wide fixed-target data 
(with p, Be, Li, and C respectively) are used to reduce the 
cold nuclear matter effects. The < pr > versus energy 
can be fitted by the function form: = p + qln^fs, 

where p, q are free parameters. The fit parameters are 
p = 0.0023 ± 0.0182, q = 0.329 ± 0.031 with x 2 /NDF = 
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41.1/15. The estimated < pt > from the fit function at 
y/s = 39 and 62.4 GeV are 1.21 ± 0.04 and 1.36 ± 0.04 
GeV/c, respectively. With these inputs, the transverse 
momentum distribution at these two cms-energies can be 
completely determined. 



FIG. 7. (color online) The ratios oi J/ip a\ | y |<i.o to a to tai as a 
function of cms-energy. The open points are the estimations 
using the two fit functions in Fig. [I] 

There are rare rapidity distribution measurements in 
p+A collisions at y/s < 200. Therefore, the universal 
energy scaling parameters of rapidity distributions are 
determined by the measurements at y/s > 200 GeV. Its 
validity at low energy (<200 GeV) range still need to 
be further investigated. But we do have various Xf dis¬ 
tribution measurements of J/ij) in fix-target experiments 
H3 m M EM EM I32H3I] • Cooperated with the a verse 
xf curve in Fig. [l] and the transverse momentum dis¬ 
tributions obtained using the strategy described in the 
above section, we can evaluate the rapidity distributions 
via the Xf distributions measurements in the fix-target 


experiments to check the validity of the rapidity inter¬ 
polation method. The ratios of J/ij) <t|i S( | <1 .o to atotai-, 
which are calculated utilizing the evaluated rapidity dis¬ 
tributions in fix-target experiments, verse cms-energy are 
shown in Fig. [7] The open points plotted in the figure 
are the estimations using the two fit functions in Fig. [4] 
In this figure, we can see that our extrapolation strategy 
also works at low cms-energy range. 


IV. SUMMARY 

We study the world-wide data of J/tp production and 
kinematics at y/s = 6.8 — 7000 GeV. We have developed 
a strategy to interpolate the J/ij) cross section, rapidity 
distribution, and transverse momentum distribution at 
any cms-energy in y/s = 6.8—7000 GeV. The rapidity and 
transverse momentum distributions measured in different 
energies have a universal energy scaling behavior. With 
this strategy, we predicted that the J/t/’ cross section 
times branching ratio at y/s = 39 and 62.4 GeV in mid¬ 
rapidity are Br(e + e~)da/dy\\y\ < i 0 = 9.04T0.69, 17.74± 
1.06 nb, respectively. 
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